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We present results from numerical solution of the Einstein field equations describing the head-on 
collision of two solitons boosted to ultra relativistic energies. We show, for the first time, that at 
sufficiently high energies the collision leads to black hole formation, consistent with hoop conjecture 
arguments. This implies that the non-linear gravitational interaction between the kinetic energy 
of the solitons causes gravitational collapse, and that arguments for black hole formation in super- 
Planck scale particle collisions are robust. 
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/. Introduction. Using ultra relativistic scattering of 
particles to probe the nature of the fundamental forces 
has a long tradition in modern physics. One reason why 
is the de Broglie relation, stating that the characteristic 
wavelength of a particle is inversely related to its momen- 
tum, and consequently probing short-range interactions 
between particles requires large momenta. In any conven- 
tional setting gravity is an irrelevant force in such inter- 
actions. However in general relativity all forms of energy, 
including momentum, gravitate, and thus at sufficiently 
high energies one would expect gravity to become impor- 
tant. The current paradigm suggests that this will hap- 
pen when center of mass energies approach the Planck 
scale, and for collisions with energies sufficiently above 
this, that black holes will be formed 

The four dimensional Planck energy Ep is ^ 10^^ GeV, 
wholly out of reach of terrestrial experiments, and as far 
as known is not reached by any astrophysical process, 
barring the Big Bang or the unobservable regions inside 
black holes. However, if there are more than 4 dimen- 
sions, intriguing scenarios have been suggested where the 
true Plank scale is very different from what is then just 
an effective 4-dimensional Planck scale [2]. A "natural" 
choice for the true Planck energy is the electroweak scale 
of ~ TeV, as this would solve the hierarchy problem. If 
this were the case, and the paradigm of black hole forma- 
tion is correct, this would imply that the Large Hadron 
Collider (LHC) will produce black holes, and that black 
holes are formed in the earth's atmosphere by cosmic 
rays Existing experimental bounds on the Planck 
energy in this context are at around 1 TeV y|. 

However, one potential problem with the above sce- 
nario, even before one considers issues regarding the 
existence of extra dimensions, physics near the Planck 
regime, etc., is whether in classical general relativity 
the generic outcome of ultra relativistic two "particle" 
scattering is a black hole for small impact parameters. 
There have so far been no solutions to the field equa- 
tions demonstrating this, and as we will outline next, 
the evidence usually presented for the case of black hole 
formation is based on a set of conjectures and the use of 
limiting-case solutions of dubious applicability. 

The main argument for black hole formation is a vari- 
ant of Thome's hoop conjecture 0: if a total amount 



of matter /energy E is compressed into a spherical region 
such that a hoop of proper circumference 2t:R completely 
encloses the matter in all directions, a black hole will form 
if the corresponding Schwarzschild radius Rs = 2GE/c^ 
is greater than i?, where G is Newton's constant and c 
the speed of light. To apply this to the collision of two 
classical spherical solitons, each with rest mass mo and 
traveling toward each other with speed v in the center 
of mass frame, let E = 27mQC^, t he total energy of the 
system, where 7 = 1/ — j (? . The largest radius R 
to be enclosed by the hoop is the rest-frame radius i^o^ 
as the Lorentz contraction only fiattens the particles in 
the direction of propagation. The hoop conjecture then 
says black holes will form if 7 ^ c^i?o/4Gmo. 

The above argument is purely classical. Quantum me- 
chanics enters with the assumption that the argument 
still holds for the collision of fundamental particles, now 
taking Rq to be the DE Broglie wavelength hc/E of the 
particle, where h is Planck's constant. Dropping con- 
stants of 0(1), the criteria for black hole formation is 
then E > {hc^ /G)^^'^ ^ the Planck energy. However, our 
goal in this paper is only to address the soundness of the 
classical hoop-conjecture argument; if it fails, there is no 
reason to expect a full quantum version to hold. 

It is not obvious that the hoop conjecture is applica- 
ble in all situations. Consider a single particle boosted 
beyond the Planck energy. Since the boosted particle's 
spacetime is a coordinate transformation of its rest frame 
geometry, there is clearly no black hole formation. As 
trivial as this example may seem, it is still insightful, as 
it illustrates that not all forms of energy gravitate in the 
same way. Here, kinetic energy, unlike rest mass energy, 
does not produce spacetime curvature, yet both forms 
of energy contribute to the mass of the spacetime, as 
measured for instance by the ADM mass \&\. The ki- 
netic energy dominates the rest mass energy by orders of 
magnitude, and for black hole formation to be a generic 
outcome the particular nature of the particles and non- 
gravitational interactions between them cannot play a 
role. Therefore, it must be the non-linear interaction be- 
tween opposing streams of gravitational kinetic energy 
that causes a black hole to form, in the process convert- 
ing kinetic to rest mass and gravitational wave energy. 

What is oft-quoted as evidence for black hole forma- 
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tion comes from the study of the cohision of two infinitely 
boosted "particles", described by the Aichelburg-Sexl 
(AS) metric [7]. The AS solution is obtained by taking a 
Schwarzschild black hole of mass m, applying a Lorentz 
boost 7, and then taking the limits 7 ^ oo and m ^ 0, 
so that the product E = 7m remains finite. The result 
is a gravitational "shock wave" , where the non-trivial ge- 
ometry is confined to a 2-dimensional plane traveling at 
the speed of light, with Minkowski spacetime on either 
side. Two such solutions, moving in opposite directions, 
can be superimposed to give the pre-collision geometry 
of the spacetime (see [H for an insightful description). 
Though the geometry is not known to the future of the 
collision, at the moment of collision a trapped surface can 
be found [9^. Assuming cosmic censorship, this would be 
an example of black hole formation in an ultra relativistic 
collision. 

There are several aspects of the infinite boost construc- 
tion that should give one pause as to its applicability to 
a large-yet-finite 7 collision of massive particles. The AS 
hmit is not asymptotically fiat, and the algebraic type of 
the metric has changed from Petrov type D (two distinct 
null eigenvectors of the Weyl tensor) to Petrov type N 
(one null eigenvector). This latter point can be thought 
of as the gravitational field changing from a Coulomb- 
like to a pure gravitational wave field. The AS metric is 
also not a good description to the geometry of a finite 
boosted particle on the shock surface; one is then left 
with the un-insightful conclusion that the description is 
good sufficiently far from the particle that its metric is 
Minkowski. It has also been argued that black hole for- 
mation is due to the strong focusing of geodesies off the 
AS Shockwave [lo|. However, there is no dynamics in this 
description (neither in the superposed AS metrics for that 
matter), and it is difficult to imagine how the geodesic 
structure can capture what is a highly non-linear and 
dynamical interaction between gravitational energy. 

To test the hypothesis that black holes form in high 
energy collisions of particles, we numerically solve the 
Einstein field equations coupled to matter that permits 
stable, self gravitating soliton solutions — these are our 
model particles. The particular soliton we use is a bo- 
son star [llj. One motivation for choosing this model 
was from earlier studies of low velocity, head-on colli- 
sions of boson stars, which suggested that as the ve- 
locity increases, gravity appears to "weaken", and the 
boson stars pass through each other exhibiting a non- 
relativistic solitonic interference pattern [l2[ . In particu- 
lar, the magnitude of the bosonic matter field developed 
interference fringes of wavelength A oc 1/P, with P the 
momentum of each boson star, which is exactly the re- 
lationship observed during the collision of Bose-Einstein 
condensates bound via Newtonian gravity. This model 
therefore seems perfect to address the genericity require- 
ment for black hole formation at super-Planck scale col- 
lisions, in that the self interaction of the matter will not 
bias the outcome toward black hole formation (as, for 
example, using black holes as model particles would). 



II. Methodology. We solve the Einstein equations, 
Rab — dab/'^R = SttTo^, usiug a variant of the generalized 
harmonic formalism [l3| with constraint damping [13] as 
described in [l5|. In this formalism, the spacetime coor- 
dinates satisfy a wave (harmonic) condition Dx^ = ^ 
with the encoding the coordinate (gauge) degrees of 
freedom. Experimentation with coordinate choices led 
us to develop what we term a damped harmonic con- 
dition, which can be written = ^ [n^ — n^] where 

is a constant, is the time-like unit normal to the 
t = const, slices, and is another time-like unit nor- 
mal field which is to be chosen so that the resulting 
coordinate system is non-singular. Such an approach 
was independently introduced in [l^, where the choice 
= {d/dt)^/a -\- \og{a / \/h)n^ ^ with a the lapse func- 
tion and h the determinant of the spatial metric, was 
proposed. We use a variant of this condition that transi- 
tions to = shortly after the collision. 

The matter is a minimally coupled complex scalar 
field (j) with mass parameter m, equation of motion 
= m?(j)^ and stress tensor Tab = 2V[a^V5]^ — 
9ab {yc<pV^4^ + m?(j)^) , where ^ is the complex conjugate 
of (j). For initial data we superimpose two boosted boson 
stars following a procedure analogous to that for binary 
black holes presented in This construction by itself 
does not fully satisfy the constraint equations. However, 
the further apart the boson stars are at t = 0, the smaller 
the error is, and we have performed tests that indicate for 
the initial separations here the error is sufficiently small 
to not qualitatively affect the conclusions. 

Each boson star is identical with a central scalar field 
amplitude 0o chose so that the maximum compactness 
2M{r)/r of each star is ^ 1/20, where the mass aspect, 
M(r), approaches the ADM mass A^ADM as r — )• 00. We 
subsequently scale all units to Mq = 2Madm- We choose 
an initial coordinate separation between the boson stars 
in the center of mass (simulation) frame of do = 250Mo, 
and give each boson star a boost of 7. Thus, with this 
compactness, the hoop-conjecture estimate of the black 
hole formation threshold is 7/^ ^ 10. 

Results. Our key result is the simple answer ^^yes!^' 
to the question of ^^do ultra relativistic boson star colli- 
sions lead to black hole formation in classical general rel- 
ativity Despite the objections we noted to the use of 
the hoop conjecture in this scenario, it therefore does ap- 
pear that the argument captures the essential physics of 
high speed soliton collisions. We find black hole forma- 
tion at 7c = 2.9 ± 10%, less than a third that predicted 
by the hoop conjecture. Of course, the latter is an order 
of magnitude estimate, and the numerical value of the 
threshold may depend on the particular soliton model. 
Note that the maximum compactness of single, stable 
boson stars is ~ 0.25 (see for e.g. [l8|), though it is not 
immediately apparent whether this is of relevance here. 

Fig. [U shows snapshots of the scalar field for several 
boost parameters at key times. For the collision be- 
ginning at rest (7 = 1), a single perturbed boson star 
forms, undergoing large oscillations that slowly damp 
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FIG. 1: Magnitude |(/)| of the scalar field from 4 different simulations, in 4 panels (left to right). The 4 sub-panels within each 
panel depict ai different times as follows (top to bottom): 1) t/Mo = 0, 2) a time at which the boson stars first completely 
overlap, 3) a short time later when |(/)| reaches a first local maximum due to gravitational focusing, 4) a late time after the 
collision. The axis of symmetry is coincident with the top edge of each sub-panel. The insets, where present, are zoom-ins 
of the central interaction regions. For the 7 = 4 case, a black hole forms near the time of sub-panel 3 — the black line in the 
corresponding inset shows the shape of the apparent horizon then, and the black semi-circle in sub-panel 4 is the excised region 
inside the black hole. 



via the emission of scalar radiation. For larger boosts, 
the initial boson star interaction exhibits the usual non- 
gravitational interference pattern, but shortly afterward 
there is some compression of the stars due to gravity. 
For the modest boost of 7 = 1.15, though the stars 
are perturbed by the compression, they pass through 
each other. Approaching the threshold with 7 = 2.75, 
the compression is much greater, and though the bo- 
son stars pass through each other the perturbation is 
strong enough to cause them to "explode". I.e., though 
the bulk of the momentum in the scalar field is concen- 
trated in two fronts propagating outward along the axis, 
a non-negligible component appears to move outward in 
spherical shells emanating from two focal points, corre- 
sponding to the locations of maximum compression seen 
in sub-panel 3. 

For the 7 = 4 case the interaction is similar to 7 = 2.75 
until apparent horizon formation; this is consistent with 
the intuition that in this regime the "matter does not 
matter" — it is the gravitational energy determining the 
dynamics, and here the scalar field is merely a tracer 
of the underlying geometry. After black hole formation, 
the resultant evolution is strikingly different. Most of the 
scalar field falls into the black hole, though a small frac- 
tion escapes. Similar to the 7 = 2.75 case, at late times 
(sub-panel 4) the matter that escapes can be traced back 
to appear to have originated as two pulses on the axis 
at the locations of maximum compression (sub-panel 3). 
However, now a piece of each outward moving wave front 
gets trapped in the black hole. The wave front remains 
connected, and with time this causes the concentric out- 



ward propagating arc patterns seen in sub-panel 4. 

Fig. [2] depicts the gravitational wave emission, as mea- 
sured by the Newman-Penrose scalar ^4, for the two 
higher 7 cases of Fig. [H Note however that the damped 
harmonic coordinates cause a strong "distortion" of the 
metric near the solitons, which eventually propagates 
outward with the gravitational wave. This prevents a 
clean interpretation of ^4, defined here with a tetrad 
aligned with the coordinate basis vectors, as representing 
the gravitational wave signal. With this caveat in mind, 
there is an interesting feature suggested by Fig. [2j The 
wavelength in the black hole formation case is consistent 
with the wave being associated with the dominant quasi- 
normal ringdown mode of the black hole. However, in the 
sub-threshold case the characteristic wavelength is quite 
a bit shorter, and there seems to be a trend of smaller 
length-scale features developing the closer to threshold 
we tune, which will continue if the threshold exhibits 
type II critical behavior [2H . As it relates to particle col- 
lisions, this suggests black hole formation may not be the 
"end of short-distance physics" , but that the low energy 
relationship equating small distances to large momenta 
ceases to be valid, and probing physics at smaller dis- 
tances requires fine-tuning the interaction energy. 
IV. Conclusions. We presented numerical results from 
a first study of the ultra relativistic collision of solitons 
within general relativity. The goal was to test if at suf- 
ficiently high energy gravity dominates the interaction, 
leading to black hole formation. We found that, for this 
class of soliton, the conjecture is true, and the threshold 
of black hole formation occurs at a boost 7c approxi- 
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FIG. 2: ^4 at t = 540Mo on the plane passing through the 
colhsion point and orthogonal to the axis, for two cases. 

mately 1/3 that predicted by Thome's hoop conjecture. 
Interestingly, a factor of ~ 1/3 also arises in calcula- 
tions of trapped-surface formation in the collision of null 
sources follo wing an S-matrix approach to the scatter- 
ing problem With 7^ = 2.9 ± 10% the ratio of 
kinetic to rest mass energy is 2 : 1; we believe this 
is sufficiently large to make a compelling case that black 
hole formation is generic in ultra relativistic particle colli- 
sions, regardless of the internal structure of the particles. 



Thus the arguments that super-Planck scale particle col- 
lisions lead to black hole formation are robust, and fur- 
thermore using black holes as the model particle to study 
the gravitational aspects of the interaction at these en- 
ergies jl, [9I, [20I is valid. Our results also suggest that 
when gravity becomes a strong player in the interaction 
at scales slightly below the Planck scale, the de Broglie 
relationship equating smaller distances to larger energies 
may cease to hold. Of course, here the nature of quantum 
gravity will be crucial, and may provide its own cut-off 
to short distance physics, though gravity will not play 
the role of the censor. 
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